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We have used the CoUider Detector at Fermilab (CDF II) to search for the flavor-changing neutral- 
current (FCNC) top quark decay t Zc using a technique employing ratios of W and Z production, 
measured in pp data corresponding to an integrated luminosity of 1.52 fb~^. The analysis uses a 
comparison of two decay chains, pp ^ tt —> WbWb —> ivbjjb and pp —f tt —f ZcWb — > iicjjb, to 
cancel systematic uncertainties in acceptance, efficiency, and luminosity. We validate the modeling 
of acceptance and efficiency for lepton identification over the multi-year dataset using another ratio 
of W and Z production, in this case the observed ratio of inclusive production ofWtoZ bosons. To 
improve the discrimination against standard model backgrounds to top quark decays, we calculate 
the top quark mass for each event with two leptons and four jets assuming it is a tt event with one 
of the top quarks decaying to Zc. For additional background discrimination we require at least one 
jet to be identified as originating from a 6-quark. No significant signal is found and we set an upper 
limit on the FCNC branching ratio Br{t Zc) using a likelihood constructed from the Ucjjb top 
quark mass distribution and the number of Ivhjjb events. Limits are set as a function of the helicity 
of the Z boson produced in the FCNC decay. For 100% longitudinally polarized Z bosons we find 
limits of 8.3% and 9.3% (95% C.L.) depending on the assumptions regarding the theoretical top 
quark pair production cross-section. 
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PACS numbers: 13.85.Rm, 12.60.Cn, 13.85.Qk, 14.65.Ha 



I. INTRODUCTION 

The standard model (SM) Lagrangian docs not con- 
tain any flavor-changing neutral current (FCNC) terms 
such as d ^ s, a consequence of its SU(2) structure 
In the SM the top quark is expected to decay via the 
charged weak current into a W boson and a bottom 
quark, t W'^b, with close to 100% branching ratio 0- 
We test this prediction by searching for FCNC interac- 
tions in top quark decays, in which the top quark de- 
cays to a Z boson and a charm quark, t — > Zc. In the 
standard model the FCNC decay t — > Zc is highly sup- 
pressed, proceeding only through radiative corrections, 
with a predicted branching ratio Br(t —^ Zc) of about 
10~^^ Q. However, some extensions of the SM (e.g. 
two-Higgs doublet models, models with extra quark sin- 
glets, technicolor models with a dynamical breakdown 
of the electroweak symmetry, etc.) predict measurable 
rates Hli. 

The production of top quark pairs, ti, is the preferred 
channel to observe the FCNC transition i ^ c at the 
Tevatron, as single top quark production has a smaller 
cross-section and much larger QCD backgrounds in the 
Zc final state. We have used data from an integrated 
luminosity of 1.52 fb~^ collected with the CDF II de- 
tector Q at the Fcrniilab Tevatron to search for events 
in which one of the top quarks decays to Zc and the 
other one decays to Wb. In order to get a sample of 
high purity, we select the leptonic decays of the Z boson, 
Z e+e" and Z IJ-~^fJ'~- In this scenario, the FCNC 
signature is most likely a pair of oppositely-charged Icp- 
tons forming a Z boson, and four jets (the b and c jets 
from the t and i, and two jets from W — > qq'), with 
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the event being kinematically consistent with the FCNC 
ti decay hypothesis. We require at least one jet with 
a displaced secondary vertex as a sign of a heavy-flavor 
quark {b or c-quark) to further suppress hadronic back- 
grounds. 

To minimize the systematic uncertainties on the par- 
ticle identification and trigger efficiencies, geometric ac- 
ceptances, and luminosity, we rely on a technique based 
on the simultaneous comparison of two decay chains: 

1. pp-^tt-^ WbWb -> ei^bjjb (sec Fig. [J), 

2. pp^ti^ ZcWb "> ilcjjb (see Fig. 

Many of the systematic uncertainties contributing to 
both decay chains are correlated and tend to can- 
cel, improving the precision and robustness of the re- 
sult. The other decay modes of top quark pairs 
(e.g._ tt WbWb liUibi2V2b, ti WbWb -> jjbjjb, 
or tt ZcWb — > £f£ic£2i^2b) have low acceptances or 
high levels of background and are not used. 




FIG. 1: A Feynman diagram for one of the processes con- 
tributing to pp tt WbWb — > Ivbjjb decay chain. 

The final states ivbjjb and ilcjjb used in this analy- 
sis contain products of the leptonic decays oi W £iy 
and Z — » ££, for which there exist precise next-to-next- 
to-leading order (NNLO) predictions of inclusive cross- 
sections multiplied by branching fractions 0]. We use a 
comparison of the measured ratio of inclusive W ^ £v 
to Z ^ ££ production to validate the lepton identifica- 
tion and trigger efficiencies in the Monte Carlo simulation 
predictions of signal and SM background to about 2%. 
This technique, which parallels that used in the search, 
will be used for precision comparisons with the standard 
model at the LHC @. 

We present a technique for measuring the background 
in the inclusive W boson sample coming from QCD 
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FIG. 2: A Feynman diagram for one of the processes con- 
tributing to pp ^ tt ZcWb iicjjb decay chain. The 
blob represents a non-SM FCNC vertex. 



multi-jet events using a data-derived model 

The 

number of mis-identified W bosons is estimated sepa- 
rately for electrons and muons by fitting the observed 
distributions in with templates from real W de- 

cays and modeled non-M^ events. 

We also present a simple technique to estimate the 
number of muons from cosmic rays in the Z boson sam- 
ple, using the distribution in the magnitude of the vector 
momentum sum, of the muon pair. This is 

an economical way of combining the usual "back-to-back" 
and momentum balance criteria for the two muons into a 
single distribution, as jX^ jJT pairs from cosmic rays have 
a very narrow peak at \P{iX^fj,')\ = GeV, while real 
Z 1^^ tJ'~ decays occupy a much larger volume in the 
3-dimensional momentum space. 

A previous limit on the branching ratio for t — > Zc Q 
is from CDF using data from Run I of the Tevatron; the 
limit is 33% at 95% confidence level (C.L.) The limit 
from precision measurements at LEP is lower, 13.7% at 
95% C.L. O. 

There is a recent CDF limit from a parallel indepen- 
dent analysis, using a different technique and a total lu- 
minosity of 1.9 fb~^, of 3.7% at 95% C.L., more restric- 
tive than the result presented here The technique 
presented here is specifically designed to reduce system- 
atic errors by using ratios of W and Z boson events, im- 
portant for much larger integrated luminosities The 
acceptance for tt— > ZcWb — > iicjjb events in this analy- 
sis is 0.303% versus the 0.43% quoted in Ref. [H. 

We derive the first upper limits on Br(t — > Zc) as a 
function of the polarization of the Z boson produced in 
a FCNC decay of a top quark. These limits cover all 
possible FCNC top-quark couplings. 

The outline of the paper is as follows. Section briefly 
describes the CDF II detector. The analysis strategy, 



which uses the SM decay modes of the top quark as well 
as the signal FCNC decay mode to allow cancellation of 
major systematic uncertainties of acceptance, efficiency, 
and luminosity, is described in Sect ion Hill Section HVl de- 
scribes the event selection, which starts with the dataset 
of events selected on central [l5l | high transverse momen- 
tum [ll| electrons and muons. The identification of jets 
containing heavy flavor is given in Section IIVDI Sec- 
tions |V] and |Vl] describe the selection of Z and W bosons, 
respectively. The modeling and validation of standard 
model vector boson production and the estimation of 
backgrounds are presented in Section fVIII Section rVIIII 
describes the technique of using the measurement of R, 
the ratio of inclusive W boson to inclusive Z boson pro- 
duction, as a check of the complex Monte Carlo samples 
generated using the detector and accelerator conditions 
accumulated over the long period of data taking. Sec- 
tion IIXI describes the modeling of the FCNC signal from 
top quark decay. The measurement of the expected SM 
contributions to the signature W boson -I- 4-jets with 
one jet identified as heavy-flavor, dominated by top quark 
pair production, is described in Section [Xj The contribu- 
tions to the reference channel, W boson -|- 4-jets, and the 
signal channel, Z boson -I- 4-jets, each with one or more 
heavy-flavor jets, are presented in Sections IXII and IXII) 
respectively. 

The estimation of systematic uncertainties on the ac- 
ceptances and backgrounds is described in Section IXIIIl 
and the limit calculations for a full range of possible lon- 
gitudinal FCNC couplings arc presented in Section IXIVI 
Section [XVI is a summary of the conclusions. 



II. THE CDF II DETECTOR 

The CDF II detector is a cylindrically symmetric spec- 
trometer designed to study pp collisions at the Fermilab 
Tevatron. The detector has been extensively described in 
the literature @. Here we briefly describe the detector 
subsystems relevant for the analysis. 

Tracking systems are used to measure the momenta 
of charged particles, and to trigger on and identify Icp- 
tons with large transverse momentu m, p x [HI . A multi- 
layer system of silicon strip detectors [3 j which identifies 
tracks in both the r — (p and r — z views [l^, and the 
central outer tracker (COT) [l3l are contained in a super- 
conducting solenoid that generates a magnetic field of 1.4 
T. The COT is a 3.1 m long open-cell drift chamber that 
makes up to 96 measurements along the track of each 
charged particle in the region \r]\ < 1. Sense wires are 
arranged in 8 alternating axial and stereo (±2°) super- 
layers with 12 wires each. For high momentum tracks, 
the COT pT resolution is CTp^/pl ~ 0.0017 GeV"^ [Hi. 

Segmented calorimeters with towers arranged in a pro- 
jective geometry, each tower consisting of an electromag- 
netic and a hadronic compartment [l9l . [20| . cover the 
central region, \r]\ < 1 (CEM/CHA), and the 'end plug' 
region, 1 < < 3.6 (PEM/PHA). In both the central 
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and end plug regions, systems with finer spatial resolu- 
tion are used to make profile measurements of electro- 
magnetic showers at shower maximum [2ll | for electron 
identification (the CES and PES systems, respectively). 
Electrons are reconstructed in the CEM with an Et 11| 
resolution of a{ET)/ET 13.5%//B^7GcV © 2% [H 
and in the PEM with an Et resolution of a{E'T)/ET — 
16.0%/\/£;T/GeV©l% [13. Jets are identified using a 
cone clustering algorithm in rj — cj) space of radius 0.4 
as a group of electromagnetic and hadronic calorime- 
ter towers; the jet energy resolution is approximately 
cr ~ 0.1 • £;t(GcV) + 1.0 GeV [H. 

Muons are identified using the central CMU, CMP, and 
CMX [23, |25] muon systems, which cover the kinematic 
region I77I < 1. The CMU system uses four layers of pla- 
nar drift chambers to detect muons with px > 1.4 GeV 
in the central region of \ri\ < 0.6. The CMP system con- 
sists of an additional four layers of planar drift chambers 
located behind 0.6 m of steel outside the magnetic return 
yoke, and detects muons with px > 2.0 GeV. The CMX 
detects muons in the region 0.6 < |?7| < 1.0 with four 
to eight layers of drift chambers, depending on the polar 
angle. 

The beam luminosity is measured using two sets of gas 
Cherenkov counters, located in the region 3.7 < \r]\ < 
4.7. The total uncertainty on the luminosity is estimated 
to be 5.9%, where 4.4% comes from the acceptance and 
operation of the luminosity monitor and 4.0% from the 
calculation of the inelastic pp cross-section [2^1 ■ 

A 3-level trigger system @ selects events for further 
analysis offiine. The first two levels of triggers consist 
of dedicated fast digital electronics analyzing a subset of 
the full detector data. The third level, applied to the full 
data from the detector for those events passing the first 
two levels, consists of a farm of computers that recon- 
struct the data and apply selection criteria for (typically) 
several hundred distinct triggers. 



mental and theoretical predictions for channels involving 
both W and Z bosons. 

We estimate R for electrons and muons separately (see 
Section rVIII[) since these particles are identified with dif- 
ferent detector subsystems. The observed numbers are 
consistent with the theoretical predictions [13, [11| and 
the previous CDF measurement [11] (see Section IVlIip . 
This cross-check was performed before measuring the 
branching ratio Br{t Zc). 

The measurement of Br{t — > Zc) is designed to be 
a measurement of the ratio between events in exclusive 
final states with a Z boson and four jets and a W boson 
and four jets. In the case of the FCNC scenario, a larger 
FCNC branching fraction leads to fewer top quark events 
decaying to -|- 4 jets, and consequently an increase 
in the rate of Z + 4 jets events. We subtract SM non- 
tt events from events with a or a Z boson and four jets 
so that the ratio is more sensitive to the FCNC signal. 
The ratio oi Z + 4 jets to -|- 4 jets increases in presence 
of FCNC events. 



IV. EVENT SELECTION 

The analysis uses events selected by the trigger system 
that contain either a central electron with Et > 18 GeV 
or a muon with px > 18 GeV [ll|. The electron dataset 
contains 75. 5M events; the muon dataset contains about 
21. 2M events. The integrated luminosity of each dataset 
is 1.52 ib-\ 

Both the observed and the simulated events (see Sec- 
tion |VlT^ are processed through the same selection cri- 
teria to identify electrons and muons, jets, W and Z 
bosons, missing transverse energy, and jets containing 
heavy flavor. Details of the selection criteria arc pro- 
vided below. 



III. INTRODUCTION TO THE ANALYSIS 
STRATEGY 

The measurement of the branching ratio of the t Zc 
decay mode is designed to be similar to the measurement 
of the i?-ratio between inclusive cross-section of W^s to 
Z's. The ratio R is defined as: 



R 



(j{W) ■ Br{W Iv) 
(j{Z) ■ Br{Z ^ U) '■ 



(1) 



where a{W) and a{Z) are cross-sections of inclusively 
produced W and Z bosons. A measurement of the R- 
ratio is itself a precise test of lepton identification ef- 
ficiencies, triggering, and Monte Carlo simulations. A 
measured i?-ratio has smaller uncertainties than a{W) 
and (j{Z) since some of the uncertainties (e.g. for inte- 
grated luminosity) completely cancel out. This makes R 
a valuable tool for precise comparisons between experi- 



A. Lepton Identification 

We use standard CDF definitions for identification 
(ID) of electrons and muons, as described below [1^. 
The same lepton ID requirements are applied to events 
from data and Monte Carlo simulations. 

The identification and triggering efficiencies for leptons 
are different for events in data and Monte Carlo, although 
they demonstrate a very similar energy dependence. To 
eliminate this inconsistency we follow the standard CDF 
practice of using correction factors ("scale factors") to 
re-weig ht the MC events (see Section ITvXSt . 

In order to maintain a high efficiency for Z bosons, 
for which we require two identified leptons, we define 
"tight" and "loose" selection criteria for both electrons 
and muons, as described below. 

To reduce backgrounds from the decays of hadrons pro- 
duced in jets, leptons are required to be "isolated" . The 
i?T de pos ited in the calorimeter towers in a cone in ?/ — 
space [la] of radius R = 0.4 around the lepton position 
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is summed, and the Et due to the lepton is subtracted. 
The remaining is required to be less than 10% of the 
lepton _Bt for electrons or px f^or muons. 

1. Electron Selection 

An electron candidate passing the tight selection must 
be central with Et > 20 GeV, and have: a) a high qual- 
ity track [131 with pt > 0.5 • Et or px > 50 GeV; b) a 
good transverse shower profile at shower maximum that 
matches the extrapolated track position; c) a lateral shar- 
ing of energy in the two calorimeter towers containing 
the electron shower consistent with that expected; and 
d) minimal leakage into the hadron calorimeter [3l|. 

Additional central electrons, classified as "loose" elec- 
trons, are required to have E^ > 12 GeV and to satisfy 
the tight central electron criteria but with a track require- 
ment of Pt > 10 GeV (rather than 0.5 ■ Et), and no re- 
quirement on a shower maximum measurement or lateral 
energy sharing between calorimeter towers. Electrons in 
the end-plug calorimeters (1.2 < \ri\ < 2.5), also classified 
as "loose" electrons, are required to have Et > 12 GeV, 
minimal leakage into the hadron calorimeter, a track con- 
taining at least 3 hits in the silicon tracking system, and 
a shower transverse shape consistent with that expected, 
with a centroid close to the extrapolated position of the 
track [H. 

2. Muon Selection 

A muon candidate passing the tight cuts must have: 
a) a well measured track in the COT [s^ with px > 
20 GeV; b) energy deposited in the calorimeter consistent 
with expectations [3J|; c) a muon "stub" (35| in both 
the CMU and CMP, or in the CMX, consistent with the 
extrapolated COT track [s^; and d) a COT track fit 
consistent with an outgoing particle from a pp collision 
and not from an incoming cosmic ray (37j . 

Additional muons, classified as "loose" , are required to 
have Pt > 12 GeV and to satisfy the same criteria as for 
tight muons but with relaxed COT track quality require- 
ments. Alternatively, for muons outside the muon sys- 
tem fiducial volume, a loose muon must satisfy the tight 
muon criteria and an additional more stringent require- 
ment on track quality, but the requirement that there be 
a matching "stub" in the muon systems is dropped. 

3. Corrections due to Modeling of Electrons and Muons in 
the MC Events 

Following the standard treatment of lepton efficiencies 
in CDF, we re- weight Monte Carlo events to take into ac- 
count the difference between the identification efficiencies 
measured in leptonic Z decays and those used in simu- 
lation [111- We then make additional corrections for the 



difference in trigger efficiencies in simulated events and 
measured in data. Corrections to trigger efficiencies are 
typically 4% for trigger electrons, 8% for trigger muons 
that traverse both the CMU and CMP systems, and 5% 
for muons in the CMX system. The average weight for 
Z e+e~ events is 0.939; for Z events it is 

0.891. 



B. Jet Identification 

Jets are reconstructed using the standard CDF cone- 
based clustering algorithm with a cone radius oi R = 0.4 
within 1 77 1 < 2.4 [39|. The jet energies are corrected 
for the //-dependent response of the calorimeters and for 
the luminosity-dependent effect of multiple-pp interac- 
tions. The simulated calorimeter response for individual 
hadrons is tuned to match that in data The raw 

energy of the jets must be greater than 8 GeV and the 
corrected energy is required to be greater than 15 GeV. 
Jets that coincide with an identified electron or photon 
are removed; i.e. each calorimeter cluster can be asso- 
ciated with either a jet, an electron, or a photon, which 
have mutually exclusive definitions to avoid any ambigu- 
ities. 

There is one case for which the jet energies are cor- 
rected to the parton level rather than to the hadron level. 
This is done to calculate the top quark mass in events 
with a Z boson and four jets (see Section IXII Ap . 

High-pT photons are not rare in hard-scattering events. 
Identifying photons as jets and then correcting them as 
jets can lead to mis-reconstructed missing transverse en- 
ergy and other kinematic variables, and can be impor- 
tant in an analysis leading to small signal samples, as 
in this analysis. Photon candidates are required to have 
no matching track with pT > 1 GeV, and at most one 
track with pT < 1 GeV, pointing at the calorimeter 
cluster; good profiles in both transverse dimensions at 
shower maximum; and minimal leakage into the hadron 
calorimeter [3l|. We require photons to be "isolated" in 
a slightly more restrictive fashion than that for the lep- 
tons: the sum of the pt of all tracks in the cone must be 
less than 2.0 GeV + 0.005 x Et- 



C. Reconstruction of Missing Transverse Energy 

Missing transverse energy (Pt) negative two- 

dimensional vector sum of Et of all identified objects in 
the event: electrons, muons, photons, jets, and unclus- 
tered energy. The unclustered energy is calculated as 
a two-dimensional vector of raw calorimeter energy cor- 
rected for the energy deposited by identified jets, elec- 
trons, muons, and photons. 
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D. Tagging of Heavy Flavor Jets 



pairs. 



We identify decays of bottom and charm quarks (heavy 
flavor, HF) with an algorithm that identifies displaced 
secondary vertexes within a jet. The primary vertex is 
identified by fitting all prompt tracks in the event to a 
vertex constrained to lie on the beam-line. Jets with 
Et > 15 GeV are checked for good quality tracks with 
hits in the COT and the silicon detector. At least two 
good tracks consistent with a common vertex arc required 
to form a secondary vertex candidate. The distance is 
calculated between the primary vertex and the secondary 
vertex candidate and projected on the jet direction. The 
jet is considered to contain a HF quark ("6-taggcd") if 
the significance of this distance is greater than 7.5a. The 
algorithm has an efficiency of approximately 50% to tag 
a 6-jet, depending on the Et of the jet, in a tt event. 
More details of the algorithm are available in [4l| . 

To model the multiple SM sources of tagged events, we 
use control samples selected from the data to estimate 
mis-tag rates (i.e. the fraction of tags coming from non- 
HF-jcts), and Monte Carlo simulated samples to get the 
contribution from SM physics processes with true heavy 
flavor jets. 

The contribution from real HF jets is estimated by 
applying the tagging algorithm to Z-t-HF and W^-fHF 
MC samples. Events with at least one fe-tag are selected. 
Each selected event is re- weighted by (1 — (1 — etag)^ 
using per- tagged-jet scale factor etag ~ 0.95±0.05 
|42| . where iVtags is the number of 6-tagged jets in the 
event, to take into account the difference in the tagging 
efficiencies between data and simulation. 

The mis-tag rate is estimated by applying the mis-tag 
parametrization (4^ to each event in a data sample that 
has all the desired characteristics except a 6-tag (called 
the "pre-tag" sample). The parametrization gives each 
jet a probability to be falsely tagged based on the jet Et, 
r], and number of tracks of good quality in the jet. 

The calculation of the mis-tag rate is performed in 
three steps. First we select all jets with Et > 10 GeV 
and \t]\ < 2.4 in the event. We then apply the mis-tag 
parameterization to the selected jets. Finally we loop 
through jets satisfying the event selection requirements 
(see Section lIVBp to calculate the probability for each 
jet to be falsely tagged as originating from a decay of a 
bottom or charm quark. The per-jet mis-tag probability 
is roughly 1%. 
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FIG. 3: The observed (points) and expected (histogram) dis- 
tributions in the invariant mass of e^e~ (upper figure) and 
fi'^ I-l' (lower figure) lepton pairs. The order of stacking in the 
histograms is the same as in their legends. 

The SM expectation for events with a Z boson and 
jets is constructed using Monte Carlo simulations of SM 
electrowcak processes such as production of WW, WZ, 
ZZ, and Z tt (see Section IVlI)) . 

The detection of Z bosons is less sensitive to the lepton 
trigger efficiencies than the detection of W bosons, since 
there are two leptons in each Z event. 



V. PRODUCTION OF Z BOSONS WITH JETS 

To be identified as a Z boson a pair of opposite-sign 
electrons or muons must have a reconstructed invariant 
mass in the mass window from 66 GeV to 116 GeV. The 
selection of Z ^ U events requires two tight leptons or 
a tight and a loose lepton. The two leptons arc required 
to be assigned the same primary vertex. Figure [3] shows 
the distributions in invariant mass for electron and muon 



VI. PRODUCTION OF W BOSONS WITH JETS 

The selection oiW ^ iiy events requires a tight central 
electron or a tight muon and greater than 25 GeV. We 
require that each W event has only one tight lepton, and 
no loose leptons. The transverse mass [4j], Mtransi^i^) , 
reconstructed from the lepton and the missing transverse 
energy is required to be greater than 20 GcV. Figure 0] 
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shows the measured and expected distributions in trans- 
verse mass for the W ^ ev and W ^lv events. 




M,,3„3(ev) [GeV] 




M„3„3(nv) [GeV] 

FIG. 4: The observed (points) and expected (histogram) dis- 
tributions in transverse mass of e + prj, (upper figure) and 
^+ Y^r^ (lower figure). The contribution from ti production is 
calculated for the case when the top quarks are decaying in 
the standard way. The order of stacking in the histograms is 
the same as in their legends. 



The SM backgrounds to events with W + jets (where 
W — > Iv) are estimated using the data and from MC 
simulations. The MC simulations are used to predict 
well-understood SM electroweak processes, such as .Z — > 
U, WW, WZ, and ZZ. Backgrounds that are largely 
instrumental, such as the misidcntification of a QCD jet 
as a Icpton from W decay, are predicted from the data. 
More details arc provided in Section IVffI 



VII. STANDARD MODEL CONTRIBUTIONS 
TO EVENTS WITH A W OR A Z BOSON AND 
JETS 

A. Monte Carlo Simulations of the Standard 
Model Processes 

The standard model expectations for the production 
of W and Z bosons are calculated from Monte Carlo 
simulations. Wc use pythia to generate W + light jets 
and Z + light jets processes and ALPGEN for generation 
of the heavy flavor processes W + HF jets and Z + HF 
jets. 

The datasets for the W and Z + light jets signatures 
are produced using a customized version of PYTHIA in 
which the px spectrum of the Z bosons, pf , has been 
tuned to CDF Run I data for < p| < 20 GeV, and 
which incorporates a tuned underlying-event [4^ and a 
requirement that Mi^^iU) > 30 GeV. The W and Z 
+ heavy flavor jets samples are produced with a ver- 
sion of ALPGEN that has built-in matching of the num- 
ber of jets from showering and matrix-element produc- 
tion [45j . Showering and hadronization of jets is done 
with PYTHIA [i^. Events from the MC generators, ALP- 
GEN and PYTHIA, are processed through the full detector 
simulation to be reconstructed and analyzed like data. 

We use the CDF version of pythia to describe the 
inclusive (i.e. before 6-tagging) production of W and 
Z bosons, used for background calculations. We use 
primarily ALPGEN samples to analyze ^-tagged events, 
as PYTHIA does not handle heavy-flavor production cor- 
rectly. ALPGEN handles radiation of additional jets bet- 
ter than PYTHIA. The inclusive production of W and Z 
bosons has only a second-order dependence on the dif- 
ference in the jet radiation of ALPGEN and pythia. ft 
has a stronger dependence on the momentum distribu- 
tion of the bosons which was tuned in the CDF version 
of PYTHIA as described above. 

The MC contributions from the SM leading order pro- 
cesses are combined into inclusive samples using weights 
proportional to the cross-sections of each contribution. 
These summed MC-samples are then compared to the 
observed events in the electron and muon decay modes 
of W and Z bosons separately. We use NNLO cross- 
sections of 2.687 nb and 251.3 pb for the production of 
W and Z bosons, respectively [231 . 



B. Electroweak Backgrounds 

Several SM processes other than DrcU-Yan production 
of W's and Z's contribute to the W and Z Icptonic sig- 
natures we use in the analysis, in particular Z — > t~^t~ , 
WW, WZ, ZZ, W Tu, and tt^ WbWb. These pro- 
cesses are estimated from corresponding MC samples, 
generated using pythia. We weight WW, WZ, and ZZ 
datasets using NLO cross-sections (13.0 pb, 3.96 pb, and 
1.56 pb, respectively [47j). 
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C. Fake Z Background from Jets Misidentified as 
Leptons 

This background consists of events in which one or 
more leptons are "fake", i.e. jets misidentified as lep- 
tons. We assume that in the samples with a vector bo- 
son and two-or-more jets, the true lepton and the fake 
lepton making up the Z in the background events have 
no charge correlation. As the number of fake Z bosons 
is small (sec below), wc use the number of same sign 
lepton pairs to estimate the QCD jet background in the 
7*/Z U sample. 

The Z — > [i^ ii~ sample, which requires 66 GeV < 
Minv(-^^) < 116 GeV, contains only 8 events with muons 
of the same sign out of 53,358 total events in the sam- 
ple. The fake muon background is consequently negligi- 
bly small. 

Same-sign electron pairs have a significant source from 
e~ pair-production by photon conversions. The ob- 
served number of same-sign electron pairs in the Z 
e+e" sample is corrected for the predicted number of 
e^e^ pairs mis-reconstructed as e'^e'^ or e"e^ using MC 
predictions for Z —>■ e+e^ production. 

Wc observe 398 same-sign electron pairs and 82,901 
e+e^ pairs. We remove the contribution of real /Z — s- 
e^e~ events from the number of observed events by sub- 
tracting the number of observed e+e~ events scaled by 
the fraction of same-sign to opposite-sign events in the 
Monte-Carlo samples for Z e+e~. The remaining 78 
same-sign electron pairs are used to estimate the QCD 
jet background in the Z e+e~ sample (see Fig. [3]). 

D. Non-W Backgrounds from Jets 

Jet production, which has a much higher cross-section 
than W or Z boson production, produces events which 
mimic the Icptonic decay of a, W boson by a mis- 
measured jet "faking" a tight isolated lepton and large 
missing energy (^x)- 

To estimate the non-W background coming from jets 
we use a data-derived model for non-W events. The num- 
ber of misidentified W bosons (non-W) is estimated sep- 
arately for electrons and muons by fitting the observed 
distributions in with templates from real W decays 
and modeled non-W events. The distributions in are 
fitted over the range < < 60 GeV using events 
that contain one tight lepton and no other leptons, with 
transverse mass Mtrans(^i') > 20 GeV (sec Fig. For 
each jet multiplicity, the non-W^ and the sum of the SM 
contributions are separately normalized in the fit to the 
distribution. The non-VF events are modeled by tak- 
ing electrons which pass all the selection criteria except 
those on the quality of the calorimeter shower (labeled 
"anti-selected-electrons" in Fig. [5]). The fractions of non- 
W events are estimated separately for events with 0, 1, 2, 
3, and >4 jets in the final state by propagating the distri- 
bution of the modeled non-W events into the region with 



> 25 GeV. The estimated fractions of non-VF events 
for each jet multiplicity are summarized in Table [J 

TABLE I: Fractions of non-W events in events with one tight 
lepton and no other leptons (inclusive W), and with >25 
GeV and Aftrans(^ + Pt) > 20 GeV. Note that this sample is 
selected without the requirement of the presence of a heavy- 
flavor jet. 

Jet Multiplicity jets 1 jet 2 jets 3 jets >4 jets 

W ^eiy + jets 0.6% 1.9% 7% 14% 20^ 
W ^ fiu + jets 0.1% 0.3% 0.9% 1.8% 2.6% 

A systematic uncertainty of 26% is assigned on the 
fractions of non-W events derived from the level of 
agreement between the shape of the data-derived non-W 
sample and the shape of '^r^ distribution of misidentified 
electrons in data. 



E. Cosmic Ray Backgrounds 

High-energy cosmic muons traverse the CDF detector 
at a significant rate and, if they intersect the beam-line, 
can be reconstructed as n^n" pairs. We remove cosmic 
ray events with an algorithm which fits the two tracks 
of the /i"*"/^" pair to a single arc composed of an incom- 
ing track segment and an outgoing segment, consistent 
in time evolution with a through-going track [33 • The 
algorithm also removes cosmic rays from events where 
only one muon is reconstructed as a W — » /i^T decay. 
It searches for hits in the COT chamber within a narrow 
road along a predicted trajectory opposite to the iden- 
tified muon. Finally, the algorithm performs a simulta- 
neous fit of the hits of the muon track and the hits in 
the predicted trajectory with a single helix to determine 
consistency with the cosmic-ray hypothesis. 

An independent estimate of the number of cosmic 
muons in the Z boson sample that have survived the 
cosmic-ray filter can be made from the distribution of 
the magnitude of the momentum vector of the 
pair, |P(/^+/i^)|. This is a simple way of combining 
the usual "back-to-back" and momentum balance crite- 
ria for the two muons into a single distribution, as cosmic 
pairs have a very narrow peak at \P{p'^ fi^)\ = 
GeV, while real Z — > fJ-^ fjL~ decays occupy only a small 
area in the 3-dimensional momentum phase space near 
P(/i+/i^)|=0. Using the |P(/i+/i^)| distribution as an 
estimator, the number of cosmic ray events in the sam- 
ple surviving the cosmic filter is negligible, as shown in 
Fig. El 



VIII. USING 7? AS A PRECISE CHECK OF THE 
MONTE CARLO SIMULATIONS 

The Monte Carlo simulation of a dataset extending 
over years, with changing detector and accelerator con- 
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FIG. 5: The measured distribution in for events with 1, 2, 3, and 4 or more jets in events with a single tight electron and 
missing energy forming a transverse mass, Mtrans, greater than 20 GeV. At least one of the jets is required to be originating 
from a heavy flavor quark. The observed distributions are compared to those from SM expectations and non-VF events (labeled 
"anti-selected-electron"- see text) in order to estimate the QCD background. The order of stacking in the histograms is the 
same as in the legend. 
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FIG. 6: The distribution |P(^+/x-)| of muon pairs from Z 
boson candidates (solid), events expected from simulations 
(stacked histogram), and from a cosmic ray sample (dashed). 
The stacked histogram is mostly Drell-Yan production; the 
other sources are considered but have negligible contribution 
to the histogram. The number of cosmic ray events in the 
search sample is negligible. 



ditions, is an exceptionally complex task, involving large 
quantities of temporal conditions stored in databases. 
Small errors in book-keeping or in properly specifying 
which data to use in the code are difficult to detect. To 
validate the modeling of the lepton identification, accep- 
tances, and triggering we use a calibration that is pre- 
dicted to better than 1.5% and is directly sensitive to 
errors affecting overall efficiencies for leptons and ^rp. 
We measure the ratio R (see Equation [T]) of inclusively 
produced W and Z bosons in their respective leptonic 
decay channels poj . 

The ratio R has been calculated at NNLO by F. A. 
Berends et al., and is predicted to be 10.67 ± 0.15 [l^. 
We measure R = 10.52 ± 0.04(stat.) using electrons and 
10.46 ± 0.05(stat.) using muons. The observed num- 
bers agree with the theoretical prediction within 2%, a 
negligible difference relative to the other systematic un- 
certainties on the t — > Zc measurement. 



multaneously two final states from decays of top quark 
pairs: pp —>■ tt —>■ ZcWb — > iicjjb and pp ^ tt 
WbWb ivbjjb, where: £ is a lepton (e or j is a jet, i/ 
is a neutrino inferred via missing transverse energy {^t^), 
b and c are "heavy-flavor" jets formed by hadronization 
of a bottom-quark or a charm-quark, respectively. This 
is done by comparing the number of expected events from 
SM ti decays and SM backgrounds to the number of 
observed events in each final state. The contributions 
from tt decays depend on two numbers: Br{t — > Zc) and 
^ti = f (pp ^ / ^'^tj where cr{pp tt) is the cross- 
section of top quark pair production at CDF and / Ldt 
is the integrated luminosity. 

Additional discrimination against SM backgrounds is 
achieved by requiring at least one of the four jets in the 
final state to be consistent with originating from a heavy- 
flavor quark (b or c quark) . The idcntiflcation of a heavy- 
flavor jet is performed with the "5-tagging" algorithm 
which is introduced earlier in Section IIVDI 

The unknown structure of the FCNC coupling is pa- 
rameterized via the polarization of the Z boson produced 
in t Zc decay, as the polarization is the only parameter 
that affects the acceptance of FCNC top quark decays. 
We vary the value of the longitudinal polarization of the 
Z bosons from 0.0 to 1.0. The final result is presented as 
a function of the longitudinal polarization. 

We reconstruct the invariant mass of the top quark, 
-^•^top, in events with two leptons and four jets assuming 
that the events are ti FCNC decays. The distribution 
of Mtop provides additional separation between standard 
model backgrounds and the FCNC signal; the top quark 
mass, Mtop, distribution for background events peaks be- 
low the FCNC signal. 



X. MEASURING TOP QUARK PAIR 
PRODUCTION IN EVENTS WITH A W BOSON 
AND FOUR JETS 



IX. THE FCNC ANALYSIS 

Assuming that the FCNC decay of the top quark t — > 
Zc is non-zero, a tt pair can decay to WbWb, WbZc, or 
ZcZc with, decay rates proportional to {1—Br{t Zc))^, 
2Br{t Zc) ■ (1 - Br{t Zc)), and Br{t Zcf , 
respectively. W and Z bosons are well identified via only 
their leptonic decay modes, which have small branching 
fractions. To keep acceptances high we require one of 
the bosons from the tt pair to decay leptonically and the 
other hadronically. 

To avoid large systematic uncertainties we analyze si- 



The measurement of the FCNC branching ratio re- 
lies on two datasets (see Section IIII[) : U + A, jets and 
I -f 4 jets, where U and i are consistent with de- 
cays of a Z boson or a boson (see Sections IVl and IVH) . 
respectively. In this section we focus only on events with 
I ^T-|-4 jets, where the majority comes from tt WbWb 
decays. At least one of the four jets in the final state is 
required to be identified as heavy flavor (HF) decay by 
the secondary vertex identification algorithm. The esti- 
mate of SM production of ly-l-HF events (e.g. W+bb) re- 
quires normalization of three key components: tt; W+bb, 
W + cc, W + c; and "non-I^" background events, which 
arise from mis-measured jet events. 
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A. Estimating the Contributions from 
it Production, Vl^+HF Production, and from non-W^ 
Backgrounds 

The dominant SM contribution to the VF+4-jet bin 
with one jet identified as heavy flavor (a "6-tag" ) is ti pro- 
duction. The production of a boson with heavy flavor, 
W+bb, W+cc, and W+c, however, dominates production 
in the W+2 jet bin. We consequently use the spectrum in 
the number of jets in W^+HF production to estimate the 
contribution from tt alone in an iterative process. We 
take the top quark pair production cross-section to be 
aitt) = 7.6 pb gi]. 

We initially assume that the fraction of non-VF events 
is negligible. We determine the normalization of the stan- 
dard model contribution to the M^-|-HF processes W+bb, 
W+c and W+cc by rescaling the respective cross-sections 
to match the total number of events observed in the W 
+ 2 jets bin. We assume that the overall normalization 
of W + bb + jets, W + cc + jets, and W + c + jets can 
be corrected by a single scale factor that is the same for 
the electron and muon channels. 

We then use this normalization of the M^+HF samples 
to estimate the remaining contribution from non-Ty's, as 
described in detail below in Section IX Bl 

We then repeat the calculation of the fraction of real 
W-I-HF events using the estimate of non-VF's, rescaling 
of the VF-hHF by a factor of 0.97 ± 0.09 to match the 
number of events in the W + 2 jets bin. The final jet 
multiplicity distributions for the W + HF sample are 
shown in Fig. [71 We find good agreement for events with 
three or more jets in the VF-I-HF sample. 

The motivation for normalizing to the two-jet multi- 
plicity bin is based on the matrix-element structure of 
associated heavy flavor production in W and Z events. 
A problem with any normalization scheme that uses the 
1-jet bin is that different diagrams contribute to the N=l 
and the N=2 jet multiplicity bins; taking into account the 
(large, particularly for charm) NLO corrections is tricky 
since the corrections differ significantly for the different 
processes. In contrast, the radiation of additional jets 
and jet matching procedures in the higher-multiplicity 
jet bins are fairly well understood (49| in comparison to 
the uncertainties in the 1-jet bin. We avoid these issues 
by normalizing the multi-jet multiplicity distribution to 
the 2-jet bin. 

We perform an additional consistency check by com- 
paring a measured top quark pair production cross- 
section with its theoretical prediction, assuming that 
there are no FCNC In the W + 4-jet bin the ratio of 
the measured cross-section to the SM expectation is 1.17 
± 0.09, where the SM background is evaluated using a 
top quark cross-section of 7.6 pb and Br{t Wb) = 
100% (i.e Br{t -> Zc) = 0). The ffr-distribution for the 
W + 4: jets events agree well with those of top quark pair 
decays (see Fig. where the contribution from the top 
quark pair production is normalized with the measured 
cross-section. The total transverse energy, Ht, is a scalar 



sum of Et of all reconstructed objects (electrons, muons, 
photons, jets, missing transverse energy, and unclustered 
energy). The transverse energies of the objects are the 
same as those used for the calculation of the missing en- 
ergy in the event. 



B. Non-W Backgrounds in the tt Sample 

The same procedure used for measuring background in 
the inclusive W bosons sample (see Section |VIID[) is used 
to measure backgrounds in the tt sample. The number of 
misidentifled W bosons (non-I^'s) is estimated by fitting 
the p-p-distribution for each jet multiplicity bin in events 
with one tight lepton and Mtrans higher than 20 GeV, 
where the transverse mass Mtrans is calculated for the 
lepton and prp. The fractions of non-M^ events obtained 
after applying the p-p-cut (^rj,>25 GeV) are presented 
in Table [Til versus the jet multiplicity (sTj i. 



TABLE II: The fractions of non-H/ QCD background (labeled 
as QCD-jets in Fig. O in events with a tight lepton (e or n), 
;Et >25 GeV, Mtrans (^ + Pt) > 20 GeV, and at least one 
6-tagged jet. 

Jet Multiplicity 1 jet 2 jets 3 jets >4 jets 
TT^ eu + jets 2.0% 4.9% 7.6% 4.7%~ 
W ^ fiu + jets 0.3% 0.9% 1.3% 2.6% 

The acceptance times efficiency, Aww^i^.^. (see Sub- 
section [XI^, is determined from the MC simulations of 
the standard model tt decays for the W + A jets bin. 
The obtained numbers are presented in Table IIIII The 
cumulative acceptances Aww-^e'^^. include the branch- 
ing fractions for W ^ Iv and W — > qq' decays. 



TABLE III: The acceptance times efficiency for I J^r^ + 4 
jets events which are produced via tt —> WbWb — > £ + 
4 jets decay chain (see Subsection IXI Aj l . The efficiencies are 
calculated from the Monte Carlo simulations and corrected to 
match the lepton identification and triggering efficiencies in 
data. 



Process 






tt^ WbWb 


e^rp-f 4 jets 


0.0128 


tt^ WbWb 


/i^Erp-f 4 jets 


0.00994 



C. Summary of the Backgrounds in W+4 Jets 



The number of events observed and the expected num- 
ber from all processes except tt production arc given in 
Tabic HVl 
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FIG. 7: The measured distributions (points) in the number 
of jets in events with a W and a 6-tag for W ^ eu and 
W — > ^v, compared to SM expectations (histogram). The 
order of stacking in the histograms is the same as in their 
legends. 



FIG. 8: The measured distribution (points) in Ht in events 
with a W and a &-tag, compared to SM expectations (his- 
togram), for the electron channel (upper figure) and muon 
channel (lower figure). The order of stacking in the his- 
tograms is the same as in their legends. 



TABLE IV; A summary of the numbers of W -|- 4 jets events. 
At least one jet in each event is required to be 6-tagged. 



Final state 


Observed 


Background (non-tt) 


ep:p + 4 jets 


252 


98.7 


/i^T + 4 jets 


219 


75.2 



XI. THE CONTRIBUTION FROM FCNC 
DECAYS OF tt PAIRS TO EVENTS WITH W/Z 
BOSONS AND JETS 

A. Acceptances for tt Decays 

We use a modified version of the madgraph Monte 
Carlo event generator [s^ to produce simulated events 
for the tt ZcWb and ti ZcZc processes, which are 
then hadronized using pythia. 

In order to calculate the rates of expected events for 



the two final states ££ + A jets and i'^r^ + 4 jets, we need 
to introduce a notation for the acceptances multiplied by 
efficiencies, {A ■ e)Y, for the decay chain "Y" of ti pairs. 
Acceptance {A ■ e)Y is a fraction of ti events observed 
in the corresponding final state. The acceptances {A ■ 
e)Y include combinatoric factors and the corresponding 
branching fractions for decays of W^s and Z's: Br{W — *■ 
iu), Br{W qq'), Br{Z «), and Br{Z -> qq). 



The acceptances {A ■ e)Y depend on the FCNC branch- 
ing ratio Br{t Zc). We divide the acceptances by 
polynomials dependent on Br{t — > Zc) to factor out the 
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terms independent of the FCNC branching ratio: 
[A ■ e)u '^ZcZc-^U+i jets 



Azw^U 
A-WZ^i^r^ 

Aww^e 
and 



Br{t Zcf 

{A ■ e)ti->ZcWb~*U+i jets 



Br{t - 


4 Zc) 


■ (1 - Br{t -> Zc)) 


(A 


■ e)ti- 


*ZcWb~*tv+i jets 


Br{t - 


^ Zc) 


• (1 - Br{t Zc)) 


{A- 




ZcWb^ef.^+i jets 



Br(t Zc) ■ (1 - Br{t Zc)) ' 

{A ■ £)tt^WbWb^£u+4: jets 



(1 -Br{t Zc))^ 



A 



(A. 



ZZ^l] 



Itt^ZcZc^l'^^+i jets 



Br{t Zc)^ 



(2) 
(3) 

(4) 
(5) 

(6) 



The values of Ay are determined using simulated samples 
where all the ti pairs decay exclusively to only one of 
the intermediate states: WbZc, ZcZc, or WbWb. The 
acceptance A-^z-tt includes two decay chains since the 
missing energy, p^, can be produced via decay W ^ tv 
or by mis-identifying Z U decay. The Z ^ U decay 
can be mistaken for a, W ^ decay when one of the 
two leptons is not identified (i.e. missing). The loss of 
the real lepton can create significant missing energy. 



B. Properties of the FCNC t Zc Coupling 

We note that the helicity structure of a possible t ^ Zc 
vertex is model-dependent. We cover the full range of 
possible helicities so as to be assumption- independent. 

The kinematic properties of i ^ Zc decay are reflected 
by the angular distributions of the decay products. This 
affects the total acceptance for the FCNC events since 
the isolation requirement is placed on all the identified 
jets and leptons. For example, the final state of the 
t — > — > ££c decay chain can be fully described by 
introducing an angle 9*, taken to be the angle between 
the direction of the top quark (anti-top-quark) and the 
positive (negative) lepton in the rest frame of the Z bo- 
son. The angular distribution of 6* has the following 
general form: 

.f{0*) = ao ■ M0*) + ai ■ /i(r) + 02 • /2(0*), (7) 

where ao, ai, and 02 are constants which depend on the 
polarization of the Z boson and whose sum is one (aq + 



J 



0.1 + ^2 = !)• The functions fi{9*) are given by: 



/o(r) = -(i-cos2(r)). 



and 



= o(i + cos(r))2 



/2(r) = -(i-cos(r))^ 



(8) 



(9) 



(10) 



The angular distribution of decay products of the 
t Wb £vb decay is parametrized with the same 
function f{9*) by taking appropriate values of the a^. In 
the case oi t ^ Wb decay the coefficients oq, ai, and 02 
are the fractions of longitudinal, left-handed, and right- 
handed helicities of the W boson, respectively. However, 
the Z boson, unlike the W boson, has both right-handed 
and left-handed couplings. Consequently, while the co- 
efficient ao is simply the fraction of the longitudinally 
polarized Z bosons, the coefficients ai and 02 are linear 
functions of the fractions of left-handed and right-handed 
helicities of the Z boson. 

The distribution of cos(0*) resulting from an arbitrary 
FCNC coupling can always be described by choosing ap- 
propriate values for the constants ai. The acceptances 
of the FCNC top quark decays Ay depend on the angu- 
lar distributions of the decay products since we require 
the isolation in a cone of 0.4 for all the identified lep- 
tons and jets. In consequence, the acceptances are func- 
tions of ao and oi (i.e. Ay = AY(ao,ai), noting that 
a2 = 1 — ao — ai). The top quark decay is symmet- 
ric with respect to the charge of the fermion {££ or qq), 
and therefore the acceptances calculated for decays of 
right-handed Z bosons and left-handed bosons are iden- 
tical. This means that the acceptances Ay can be fully 
parametrized with the fraction of longitudinally polarized 
Z bosons (i.e. Ay = Ay{ao, 1 — ao) = AY(ao))- 

We compute each acceptance Ay for five values of 
the fraction of longitudinally polarized Z bosons using 
Monte Carlo simulated events. This allows us to calcu- 
late the acceptances Ay for any fraction ao by interpo- 
lating the acceptances Ay between the points measured. 
The acceptance Aww^i'?,^, is a constant since it does 
not have any FCNC vertices. The other acceptances, 
Azz^ee, Azw^ee, Awz^e^^T,, and Azz^e^^, have hncar 
or quadratic dependences on the fraction of the longitu- 
dinal helicity of the Z bosons: 



A 



zz- 



(ao) = al • + 2 . ao • (1 - ao) • ATd-.u + (1 - ^o)^ ' A 



Azw^u(ao) — 0-0 ■ A^zw~*ee + (1 ^ 'lo) ' A^zw- 



left 

zz^u^ 



(11) 

(12) 
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Awz~,ef.T,{ao) = ao ■ A\^§^^^^ + (1 - ao) • ^wz-^^g^, (13) 

and 

Azz^i^^iao) = al ■ A^z'z-^e^^ + 2 • ao • (1 - ao) • ^zT-^^ + (1 - a-of ■ A^z%^i^^, (14) 



where Ay"^ are measured for the longitudinally-polarized 
component of the Z decays, ^y^' are for the left-handed 
component, and the value of A™" is obtained using 
FCNC events where the Z bosons are mixed with 50% 
left-handed and 50% longitudinal polarizations. The ac- 
ceptance Azz^u has a quadratic dependence on ao since 
it accounts for the two FCNC decays of the top and anti- 
top quarks. The numerical values of the acceptances are 
tabulated in Sections |X] and IXIIl 



XII. MEASURING THE CONTRIBUTIONS 
FROM FCNC AND SM PROCESSES IN EVENTS 
WITH A Z BOSON AND FOUR JETS 

At this stage we consider only events which have two 
leptons consistent with a parent Z boson and at least 
one 6-tagged jet. We use the jet multiphcity distribu- 
tion (see Fig. [9|) to constrain the number of non-SM 
Z+^-]et events. We do this by scaling the total Z+W 
component, ZA-q {q^c,b), to the number of (observed - 
mis-tagged) Z+2 jets events in the electron and muon 
modes simultaneously. The fraction of mis-tagged, Z+q 
(q = u,d,s), events is estimated from data using inclu- 
sive Z + jets events (see Section llVDp . The number of Z 
+ 4 jets events observed and the expected number from 
all SM processes are given in Table fVl 



TABLE V: A summary of the numbers of Z + 4 jets events. 
At least one jet in each event is required to be 6-tagged. 
Final state Observed SM Background 



+ 4 jets 



8.4 



fi + 4 jets 



6.9 



The FCNC signal contribution is divided into two 
parts, ZcWb and ZcZc, since the 6-tagging rates are dif- 
ferent. We summarize the acceptance and the efficiency 
measurements for the Z + 4 jets channel in Tables IVIl 
IVIIl IVIIIl and IIXI The case when a leptonic decay of 
a Z boson is misidentified as the leptonic decay of a 
boson is taken into account in Tables IVIIIl and IIXI 

The top quark mass A/top is used as a discriminating 
variable against the SM backgrounds as was mentioned 
earlier in Section HXl More details on the calculation of 
Mtop are provided later in Section IXII Al We recalculate 



the acceptances ^^z^w ^^'^ ^^zw^ee ^'^^ '^^ 
the top quark mass distribution by multiplying a cumu- 
lative acceptance Ay (Y is ZZ U or ZW — > U) and 
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FIG. 9: The measured distribution (points) in the number of 
jets in events with a Z and a fe-tag, compared to SM expecta- 
tions (histogram) , for the electron channel (upper figure) and 
muon channel (lower figure). We normalize to the average 
of the Z — > e^e~ and Z 2 jet bins. The order of 

stacking in the histograms is the same as in their legends. 



the fraction of events in the i'th bin: 



A- 



Y 



A 



Y 



E 

k 



Nk 



(15) 



The obtained acceptances Ay depend on the recon- 
structed top quark mass of the ti FCNC events. 
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TABLE VI: The acceptance times efficiency for the dilepton 
signature from the inclusive FCNC decay of ti^ ZcZc 
ii + ccjj for different values of the longitudinal fraction of Z 
bosons, for electron pairs and muon pairs separately. The SM 
branching ratios for the Z —t U decays are included. 

Process Azz-,u 

Longitudinal fraction is ao=0.00 
tt^ ZcZc -> e+e~ + 4 jets 0.00185~ 
a-> ZcZc + 4 jets 0.00178" 

Longitudinal fraction is ao=0.50 
tt^ ZcZc e+e~ + 4 jets 0.00203~ 
tt~* ZcZc ~> ^1+^- + 4 jets 0.00192" 

Longitudinal fraction is ao = 1.00 
tt^ ZcZc ^ e+e~ + 4 jets 0.00222" 
tt^ ZcZc ^1+^- + 4 jets 0.00205" 



TABLE VIII: A summary of the acceptance times efficiency 
for the contribution to the single lepton+p,j, signature from 
the inclusive FCNC decays of tt^ WbZc ^ + bcjj (i.e. 
the decay of a Z boson is mis-identified as the decay of a 
W boson) and tt— > WbZc —> £v + bcjj. Standard model 
branching ratios are included. The acceptance Awz^i^^, is 
the sum of acceptances for the decay modes which contribute 
to the signature of ^ + + 4 jets. 



Process 


Awz-*i 


Longitudinal fraction is 


ao^O.OO 


tt^ WbZc eu + 4 jets 


0.00927 


tt^ WbZc ^ + 4 jets 


0.00179 


tt^ WbZc ^ ^i^ + 4 jets 


0.007915 


tt~* WbZc ^ /.tpx + 4 jots 


0.002180 


Longitudinal fraction is 


ao^l.OO 


tt^ WbZc eu + 4 jets 


0.00967 


tt^ WbZc ^ +4 jets 


0.00185 


tt^ WbZc -» ^ii/ + 4 jets 


0.00817 


tt-* WbZc ^ fi^rj, + 4 jets 


0.00227 



TABLE VII: A summary of the acceptance times efficiency 
for the dilepton signature from inclusive FCNC decays of 
tt— > ZcWb II -\- bcjj, for different values of the longitu- 
dinal fraction of Z bosons, for electron pairs and muon pairs 
separately. 



Process 


Azw^u 


Longitudinal fraction is ao 


=0.00 


ZcWb e+e" + 4 jets 


0.00275 


ZcWb 11+ fi' + 4 jets 


0.00267 


Longitudinal fraction is ao 


=1.00 


tt^ ZcWb ^ e+e" + 4 jets 


0.00313 


ZcWb 11+ fi' + 4 jets 


0.00293 



A. Fitting the Top Quark Mass 



TABLE IX: The acceptance times efficiency for the contri- 
bution to the single lepton-|-px signature from the inclusive 
FCNC decay of tt— > ZcZc —> i + + ccjj, where at least 
one di-leptonic decay of Z boson has been mis-identified as 
the decay of a boson. Standard model branching ratios 
are included. Note that this channel depends on the square 
of the FCNC branching ratio for the Z, and so its contribu- 
tion is suppressed relative to that from the case where only 
one Z decays by FCNC (see Table rvm)l . 



Process 




Longitudinal fraction is ao 


=0.00 


tt—^ ZcZc ^ e + -1- 4 jets 


0.000873 


it^ ZcZc ^+ + 4 jets 


0.00127 


Longitudinal fraction is ao 


=0.50 


tt— > ZcZc e+ + 4 jets 


0.000858 


tt-> ZcZc ^1+ ^x + 4 jets 


0.00132 


Longitudinal fraction is ao 


=1.00 


tt—^ ZcZc ^ e + ^x -1- 4 jets 


0.000838 


tt^ ZcZc ^ ^ + ^x + 4 jets 


0.00137 



We reconstruct the value of Aftop for each candidate 
event that contains at least two leptons consistent with 
a parent Z boson and at least four jets. The procedure 
is very similar to that of the CDF top quark mass mea- 
surement I53l. 



The value of Mtop is calculated by minimizing the 
distribution, which is based on the assumption that the 
event is pp tt ^ Z + 4 jets ££ + 4 jets. The 
minimization takes into account every combination of 
the jets in the event since we do not know the true 
jet-parton assignments. To do so we loop through all 
possible permutations and select the one with the low- 
est x^- The top quark mass distribution obtained for 
tt ZcZc ~* U + 4 jets decays does not differ signifi- 
cantly from that of WbZc decay. The exact formula for 



the has the following structure: 

X^(M,op)= E ^^^^^ 



f 1,^2 Jots ' 

"i ^ W 



E 

x,y 

{M{l+l-)-Mzf , {M(W +j)-M,o^f 



^ z 



^ top 



(A/(Z+j)-Mtop)^ 

^ top 



(16) 



The first term contains the fitted transverse energies of 
the leptons and four jets within the corresponding exper- 
imental resolutions. The second term includes the x- and 
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y- components of the unclustered energy. The expression 
also contains terms for the reconstructed masses of the 
W , Z, and the two top quarks (i.e. t Zc iijet 
and t — > Wb — » 3 jets). The function inchides all 
the top-specific corrections of jet energy scales and en- 
ergy resolutions used in the singie-lepton top quark mass 
measurement [53j . 

We process the Z + A jets events from data and sim- 
ulation samples with the same top quark mass fit com- 
puter code so that we can compare the Mtop distributions 
between data, the SM expectations, and a hypothetical 
FCNC signal. The comparison is shown in Fig. [TOl the 
data agree well with the SM background distribution. In 
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FIG. 10: The measured distribution (points) in the fitted top 
quark mass in events with a Z and four jets with at least one 
b-tagged jet, compared to the SM expectations and an FCNC 
signal (stacked histogram), for the electron channel (upper 
figure) and muon channel (lower figure) . The branching frac- 
tion for the FCNC signal is taken from Table IXlIl The order 
of stacking in the histograms is the same as in their legends. 

the following section we describe the evaluation of the 
systematic uncertainties that go into making this state- 
ment quantitative and setting a limit on a FCNC signal. 



XIII. SYSTEMATIC UNCERTAINTIES 

We discuss separately the systematic uncertainties in- 
volving the acceptances and backgrounds in the following 
two subsections. 



A. Systematic Uncertainties on the Acceptances 

The uncertainties on the five acceptances used in de- 
termining the limit. Ay, defined in Section [IXI are sum- 
marized in Table [Xj For each of the Ay the effect of un- 
certainties in the jet energy scale, initial and final state 
radiation, lepton identification efficiencies, parton distri- 
bution fimctions, and the identification ("tagging") of 
bottom quarks and charm quarks have been taken into 
account. 

To estimate the systematic uncertainty from each of 
these sources, we vary each of the parameters listed be- 
low by one standard deviation (icr) and recalculate the 
acceptances Ay- The effect of the uncertainty for each 
of the sources is correlated among the Ay, and these 
correlations are taken into account in the limit-setting 
procedure. 

Of the two acceptances that contribute to the sig- 
natures containing two charged leptons, Azz^u and 
Azw^ee, the latter dominates as it depends linearly on 
the FCNC branching ratio of the Z boson, while the con- 
tribution corresponding to Azz-*u enters as the square. 
In a similar fashion, the single lepton -I- signature is 
dominated by the SM decay of the top quark pair into 
W^W~bb, with an acceptance Aww^e^^,, as there is 
no FCNC branching ratio in the rate. The process de- 
scribed by Awz^i^rj. is suppressed by a single factor 
of the FCNC branching ratio, while that described by 
Azz—te^,j, is quadratic, and hence makes a very small 
contribution. 

The largest systematic uncertainties in the dominant 
processes in the dilepton and singie-lepton modes are the 
uncertainties in the efficiency for identifying b- and c- 
quarks. For 6-quarks, we follow the prescription used 
previously in CDF studies of the top quark, and use a sys- 



tematic uncertainty on the tagging efficiency of 5% 42 1. 
Similarly, for c-quarks, we assign a 15% uncertainty 41 1. 

The next largest contribution to the systematic un- 
certainties is from uncertainties in the calibration of jet 
energies (53j . The systematic uncertainties are positively 
correlated for all the Ay. 

The contributions from lepton identification and trig- 
ger efficiencies are limited by the precision check of the 
_R-ratio (see Section IVIIip . Wc assume that the recon- 
struction and the triggering efficiencies of electrons and 
muons are not correlated. We note that acceptances and 
trigger efficiencies are correlated for W ^ iiy and Z —> ££ 
decays to leptons of the same flavor. This means that 
Azz—fU would be mis-estimated by the same percentage 
as Azw-tie for leptons of the same flavor. The same 
holds true for Awz~*e^^, Aww^e^.^, and Azz^if,^.- 
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TABLE X: A summary table of the systematic uncertainties on the acceptances. Correlations are taken into account in the 
calculation of the limit. The abbreviation "lept." stands for the systematic uncertainty due to lepton identification and 
triggering. 



Systematic Uncertainty in % 




^i^zw^ee) 


''{^WZ^i ) 






'^ZZ^U 










JES 


2.5 


2.6 


2.7 


2.4 


6.4 


ISR 


0.5 


0.5 


0.5 


0.5 


0.5 


FSR 


0.6 


0.6 


0.6 


0.6 


0.6 


PDFs 


0.9 


0.9 


0.9 


0.9 


0.9 


HF quark ID 


10.2 


5.0 


5.0 


4.1 


10.6 


ID and triggering of electrons 


1.6 


1.6 


0.1 


0.1 


0.1 


or 


ID and triggering of muons 


2.8 


2.8 


0.9 


0.9 


0.9 


Total 


10.6 © lept. 


5.8 © lept. 


5.8 © lept. 


4.9 © lept. 


12.4 © lept. 



The systematic uncertainties in the Ay due to lepton 
identification and triggering arc estimated using devia- 
tions between the measured cross-sections of inclusive 
Ws and Z's, used in calculating the ratio i?, from their 
theoretical values: 



Acr(Z ^ U) 

a{Z ££) 



Azz^i 



5{Azw-*u) 
Azw^u 



(17) 



and 



5{A 



WZ^i 



i{Aww-*i'^^. 

Aww^i^j. 
5{Azz^e^^) 



A 



zz- 



(18) 



The uncertainty on the integrated luminosity does not 
contribute at the first order to the measurement of 
Br{t — > Zc) since it is positively correlated between 
a{W iv) and a{Z U). 

The deviation of the measured R ratio is 



a{Z 
Ao-(W ^ 



U) 



£v) 



\a{Z U) 
Acr{Z U) 



(t{W £v) a[Z U) 



(19) 



Therefore, the connection between the deviation in the 
i?-ratio and the uncertainties of the Ay is the following: 



Ai? 



Azw^u 



I 



A 



zw^u 



^Aww-^e^rj. J \Aww-*ef 
S{Azw->u) _ S(Aww^e^^) 
Azw^u 



(20) 



We treat 



S{A 



and 



as negatively cor- 



Aww^a^rj. Azw- 

related as it is the most conservative case. Also this 
treatment insures the constraint from the i?-ratio. 

Contributions from other sources are significantly 
smaller than those from heavy flavor identification and 



jet-energy scale. The effect of initial and final state radi- 
ation (ISR and FSR) on Aww-^f^^ was studied in Ref. 
[5^ . We expect that FSR will contribute to the uncer- 
tainties in the other three Ay in the same way since we 
require four jets in the final state for all four channels 
and the samples are triggered on leptons. The ISR un- 
certainty should also contribute identically to the uncer- 
tainties of the four acceptances Ay- The uncertainties 
are found to be 0.5% for ISR and 0.6% for FSR, assumed 
to be 100% correlated across all Ay- 

The uncertainties arising from parton distribution 
functions (PDF's) can also propagate into the accep- 
tances. However, the dominant effect of changes in the 
PDF's is on the production of the tt pairs and not on 
the decay kinematics. The effect of the uncertainties was 
also studied in Ref. [5J|. The total uncertainty is 0.9% 



and is 100% correlated for the four Ay- 



B. Systematic Uncertainties of the Backgrounds 

The sensitivity of this search for a Z boson and a charm 
quark coming from top quark decay depends strongly on 
the understanding of SM W boson and Z boson produc- 
tion in conjunction with heavy flavor {W /Z+YiY). We 
summarize the systematic uncertainties of backgrounds 
in both the single lepton and di-lcpton signatures (the 
terms Bii, and Bgg in Equations and in Table IxH 
and discuss them below. 

The largest uncertainty in the background comes from 
modeling the production of W bosons and Z bosons 
accompanied by heavy-fiavor and additional jets. The 
Z+HF and W^-I-HF backgrounds are modeled by ALP- 
GEN [H] , and hadronizcd with pythia [111 . The predic- 
tions suffer from uncertainties in the modeling procedure. 
In particular, the expected number of events in the W / Z 
+ 4 jets category enters directly into the calculation for 
the final result. To make an estimate of the uncertainty 
on the expected number of W/Z+E.F events, we assume 
that there is a set or parameters which allows ALPGEN 
to model the data precisely. A deviation from the "ideal 
set" can be estimated using inclusive Z + jets events 
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TABLE XI: The relative systematic uncertainties (%) on the 
backgrounds for 4-jet semi-leptonic and dilepton final states 
of tt pairs. The contributions from uncertainties in the Monte 
Carlo modeling and in the rate of misidentified heavy-flavor 
jets (mis-tags) are (conservatively) taken to be correlated in 
the computation of the limit. 



Systematic Uncertainty in 


% 


+ 4jets 


M + 4jets 


W/Z+HF+3ets 




20 


20 


MC Modehng 








Mis-tags 




15 


15 


W/Z+RF 




2.5 


8 


Normalization 









with jet multiplicity below three. A comparison between 
data and ALPGEN simulations is shown in Fig. [TlJ The 
observed deviation on the rate of radiation of one extra 
jet in the inclusive sample is less than 5%. We assume 
independent gluon emission, and so take 10% as the es- 
timate of the uncertainty on this ALPGEN prediction for 
the radiation of 2 extra jets in the inclusive sample. How- 
ever, the slopes of the N-jet distribution arc predicted to 
be different in the inclusive and HF samples, with the 
factors for each additional jet being 5.0 and 2.7 in the 
inclusive and 5-tagged samples, respectively. The ratio 
of 5.0 to 2.7 makes a relative difference of 1.85 between 
radiating an extra jet in inclusive and tagged samples. 
Wc consequently increase the 10% deviation by a factor 
of 2 (rounding 1.85 up), to 20%. 

The sensitivity of the limit to the number of 4-jet 
Z+HF events was calculated by performing a set of 
"pseudo-experiments" with different levels of the system- 
atic uncertainties of the backgrounds. The limit was 
calculated with this systematic uncertainty set to zero, 
set to 20% (nominal), and set to 40%. The respective 
shifts in the limit are -0.1%, zero (by construction), and 
-fO.1%, respectively. The weak dependence is caused by 
the measurement technique; we measure a ratio of top 
quark events between Z + A jets and W + 4 jets final 
states. An increase in the number of background events 
leads to a decrease in the ti cross-section measured with 
W + A jets events. Simultaneously it leads to a decrease 
in the upper limit on the number of FCNC signal events 
in the Z + A jets final state. 

The method of predicting misidentified heavy flavor 
(mis-tags) by applying a parameterization of the rate for 
a light-quark jet or gluon jet being mis-idcntificd as a 
jet from a charm or bottom quark to jets in a sample 
before heavy flavor identification contributes a significant 
systematic uncertainty to the background estimates. We 
vary the mis-tag probability calculated by the standard 
CDF algorithm used in the measurement of the top quark 
cross-section [53] jet-by-jet by ±15 % (i.e. a factor of 0.85 
or 1.15) to estimate the contribution to the uncertainty. 

A smaller contribution to the uncertainty is due to 
the overall normalization of the predicted SM boson-|-HF 
contribution. The normalizations of the background dis- 
tributions from VK+HF and Z+HF events arc treated 
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FIG. 11: The measured distribution (points) in the number of 
jets in events with an inclusive decay of ^ — > jj,'^ ii~ , compared 
to SM expectations (stacked histogram). The order of stack- 
ing in the histogram is the same as in the legend. The Z-|-jets 
processes (Drell-Yan, Z+b, and Z+c) are modeled with ALP- 
GEN. The lower plot shows the difference between the data 
and predictions. 



as independent, and are chosen to match the number of 
observed events in the VF-l-HF-l-2 jets and Z-|-HF-f 2 jets 
channels, respectively, as discussed in detail in Sections 1X1 
and IXIII The finite statistics of the 2-jet bin of the data 
contributes an uncertainty of 2.5% to the single lepton 
and di-lepton signatures, respectively. 



The 6% uncertainty of the measured luminosity affects 
only processes that are normalized absolutely: WW, 
WZ, and ZZ production. Consequently, the contribu- 
tion from the uncertainty of luminosity to the final result 
is negligible (< 0.1%). 
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XIV. STATISTICAL EVALUATION OF THE 
LIMITS ON Br{t -» Zc) 

At this point we have all the ingredients needed to 
evaluate limits on the FCNC branching ratio Br{t 
Zc). The branching ratio is evaluated by comparing the 
numbers of expected and observed events in two final 
states, U+ijcts and £ ^rp+4jets, using Poisson statistics. 
The numbers of observed events are denoted as Ni^ and 
N(( for final states £ ^x+^jets and ^^+4jets, respectively, 
the numbers of expected events are denoted as Xi^, and 

To avoid large systematic uncertainties we simultane- 
ously analyze two final states from decays of top quark 
pairs: pp ^ tt ^ ZcWb — > ££cjjh, and pp tt ^ 
WbWb £ ^rj,bjjb. This is done by comparing the 
number of expected events from SM tt decays and SM 
backgrounds to the number of observed events in each 
final state. The contributions from tt decays depend on 
two numbers: Br{t Zc) and Ntt = cr{pp — > ti) J Ldt, 
where (7{pp — tt) is the cross-section of top quark pair 
production and / Ldt is the integrated luminosity. We 
treat Br{t — > Zc) and iVjf as free parameters in the cal- 
culation of the limit on the FCNC branching ratio. The 
result of the comparison is presented as a likelihood which 
is a two-dimensional function of Br(t — > Zc) and Ntt- 
We use the likelihood distribution to estimate limits on 
the FCNC branching ratio Br{t Zc) using a Bayesian 
approach. 

For simplicity, let us consider the case in which we 
observe only two categories of events: Ni^ and Nu, by 
applying some set of selection requirements. Later we 
will show how to generalize this approach to be used with 
more categories of selected events. This is done since we 
will consider events with electrons and muons separately 
and we use a binned distribution of Mtop for ££+A jets 
events. 

We assume that the top quark has only the two decay 
channels Wb and Zc, and so Br{t Wb) + Br{t 
Zc) = 1 . The number of expected ti pairs is 

Ntt = o{vV ^tt)- J Ldt, (21) 

where cr^pp — > tt) can be taken a priori since it is inde- 
pendent of any FCNC physics. 

The expected numbers of events in each of the decay 
modes are estimated as follows, where we use the nota- 
tion Bz = Br{t Zc): 

Xi^ w Bii, + NtiAww^i'^^ (22) 

and 

Xu ~ Bu + NtiAzw~^u ■ Bz- (23) 

The complete formulas are presented in [57| . In the for- 
mulas above Bi^, and Bu are non-top SM contributions 
(backgrounds) to final states £ ^-p -I- 4jets and ££ + 4jets, 



respectively; Ay is acceptance for a decay mode "Y" (see 
Section HX]). 

The limit on the ratio Br{t — > Zc) is estimated using 
probability density (i.e. likelihood) function defined as: 

L{Bz, Ntt) = P{Ni,, Nu\Bz,Ntt) (24) 

i.e. 

L{Bz, Ntt) = P{Ne,\Xe,)P{Nu\Xu), (25) 

where 

X'^e-^ 

P{N\X) = (26) 

is a Poisson distribution. The likelihood L{Bz,Ntt) is 
defined in the physical region of parameters Ntt 1^ and 
< < 1. 

The complete set of systematic uncertainties is in- 
cluded in the likelihood function using a Monte Carlo 
simulation which takes into account the correlations be- 
tween the uncertainties. 

To discriminate the FCNC signal from the expected 
SM background, we use the distribution in the recon- 
structed top-quark mass, Mtop, for Z+4 jets events. 
Events from the signal process should form a distinguish- 
able peak at the top quark mass. We combine probabil- 
ities for each bin of the reconstructed top-quark mass 
distribution 

\{P{Nl,\Xl,), (27) 

i 

where the index i refers to the i'th bin of the distribu- 
tion in the top mass. This requires calculating the ac- 
ceptances A^z^^^ and A^^-^^^^ for each bin of the re- 
constructed top quark mass histogram. 

We note that the electron and muon decay modes of 
the top quarks are treated separately up to this point of 
the analysis in order to better understand the systematics 
of both. The two channels are then included together in 
the final likelihood function L{Bz, Ntt)- 

The likelihood hmction is used to construct a posterior 
probability density P{Bz\DATA), where DATA refers to 
the numbers of observed events, Nii, and N'^^, in the elec- 
tron and muon channels {£ = e oi £ ^ fi). The posterior 
probability density function is converted into a limit on 
the FCNC branching ratio Br{t — > Zc) using a Bayesian 
approach. 

A. Numerical Computation of the Likelihood 
Distribution Function L(Bz,Ntt) 

The observed distribution of the likelihood (computed 
for t ^ Zc decays where the Z bosons are 100% longitu- 
dinally polarized) is presented in Fig. [T^l 

A likelihood distribution is calculated for each given 
value of helicity of the t ^ Zc coupling since the ac- 
ceptances Ay vary for different structures of the FCNC 
coupling. 
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The distribution for P{Bz\T)ATA), calculated for 
100% longitudinally polarized Z bosons, is shown in 
Fig.Ill 
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FIG. 12: The likelihood distribution L{Bz,Ntt) calculated 
as a function of Nft and Bz ~ Br{t Zc). The distribu- 
tion is for FCNC decays oi t ^ Zc with 100% longitudinally 
polarized Z bosons. 



B. Computation of the Posterior 

P{Br(t Zc)|DATA) 
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FIG. 13: The distribution for P{Br{t Zc)|DATA), calcu- 
lated for 100% longitudinally polarized Z bosons using Guas- 
sian prior. 



The posterior probability density functions 
P(-B2 1 DATA) are computed from the likelihood 
bmctions L{Bz,Nii) using a Bayesian approach as 
follows: 

P(DATA|Bz)= / L{Bz,Nu)-MNu)dNu (28) 

JO 

P(i^.|DATA)= m^^MBz):.iBz) ^^^^ 

[ PiDATA\Bz)-T:iiBz)dBz 
Jo 

where 7ro(A^tt) is the a priori probability density function 
of iVjf and tti (Bz) is the a priori distribution of Bz which 
is taken to be flat in the physical region (it is 1.0 for 
< -Bz < 1 and zero everywhere else). The distribution 
of Tro{Nft) represents the prior knowledge of the top pair 
production cross-section, a{pp — > tt). 

We consider two choices of the iroiN^t) prior distribu- 
tion: flat and Gaussian. The flat distribution does not 
contain any information regarding the theoretical predic- 
tions of cr{pp —>■ tt). It is just a constant. The Gaussian 
distribution is derived using the theoretical estimates of 
the top pair production cross-section cr[pp tt) [sst and 
the integrated luminosity. The theoretical estimate of the 
top pair production cross-section is presented as a func- 
tion of top quark mass Mtop. The measured top-quark 
mass is 170.9 ± 1.8 GeV The luminosity is 1.52 ib-\ 
with an uncertainty of 6%. The Gaussian prior allows us 
to take into account the theoretical FCNC-independent 
knowledge of a{pp — > tt). 



C. Computation of the Upper Limits on Br{t — > Zc) 

We use the posterior function P(i?z [DATA) to calcu- 
late the upper limit i?^™ on Br{t Zc) (i.e. Bz) by 
solving the equation: 

P = I P{Bz\BATA)dBz, (30) 
Jo 

where /3 is 0.95 (95% C.L). The upper limits versus the 
helicity of the Z boson are summarized in Table IXIIl 

We perform statistical cross-checks of the measured 
upper limits using pseudo-experiments. The pseudo- 
experiments are generated randomly assuming that there 
is no contribution from FCNC processes, i.e. by setting 
Br{t Zc) = 0. The expected upper limit for 100% 
longitudinally polarized ^'s on Br{t Zc) is 8.7±2.6%, 
consistent with the observed limit of 8.3%. 



XV. CONCLUSIONS AND RESULTS 

Taking into account systematic uncertainties on Monte 
Carlo simulations, 6-tagging, mis-tag modeling, and lep- 
ton identification, we find an upper limit at 95% C.L. on 
the branching ratio of t — > Zc of 8.3% for FCNC decays 
where the Z bosons are 100% longitudinally polarized. 
The result is primarily statistics-limited. It can be sig- 
nificantly improved with more data if the number of Z+A 
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jets events is high enough to do a shape analysis of the 
top quark mass distribution. 

To be assumption-independent we parametrize the 
Umit on Br{t — » Zc) as a function of the fraction of 
longitudinaUy polarized Z bosons. The parametrization 
allows us to cover the full range of all possible helicity 
structures of the t — > Zc vertex. The upper limits arc 
calculated at 95% C.L. for five fractions of longitudinally 
polarized Z's using 1.52 fb~^ of data. The results are 
presented in Table IXIII for both the Gaussian and the 
flat priors. The limits vary between 8.3 and 9.0% for the 
Gaussian prior depending on the polarization of the Z 
boson and are about 1% less restrictive for the flat prior. 



TABLE XII: The upper limits on the FCNC branching ratio 
Br{t — » Zc) in % as a function of the longitudinal fraction 
of the Z bosons in the FCNC coupling {t -> Zc) at 95% CL. 
The limits labeled Gaussian prior use as input the theoretical 
cross-section of u{pp — > tt); the limits labeled Flat prior are 
theory-independent . 



Longitudinal Fraction 0.00 


0.25 


0.50 


0.75 


1.0 


Gaussian prior 9.0% 


8.8% 


8.6% 


8.5% 


8.3% 


Flat prior 10.2% 


10.0% 


9.7% 


9.5% 


9.2% 



XVI. ACKNOWLEDGMENTS 

We thank the Fermilab staff and the technical staffs of 
the participating institutions for their vital contributions. 
We thank Mary Heintz for unfailing computer support. 
We are grateful to Michel Herquet for explaining to us 
how to incorporate the FCNC couplings into madgraph 
code, and then installing the model in madgraph. 

This work was supported by the U.S. Department of 
Energy and National Science Foundation; the Italian Isti- 
tuto Nazionalc di Fisica Nucleare; the Ministry of Educa- 
tion, Culture, Sports, Science and Technology of Japan; 
the Natural Sciences and Engineering Research Council 
of Canada; the National Science Council of the Repub- 
lic of China; the Swiss National Science Foundation; the 
A. P. Sloan Foundation; the Bundesministerium fiir Bil- 
dung und Forschung, Germany; the Korean Science and 
Engineering Foundation and the Korean Research Foun- 
dation; the Science and Technology Facilities Council and 
the Royal Society, UK; the Institut National dc Physique 
Nucleaire et Physique des Particules/CNRS; the Russian 
Foundation for Basic Research; the Ministerio de Cien- 
cia e Innovacion, and Programa Consolider-Ingenio 2010, 
Spain; the Slovak R&D Agency; and the Academy of Fin- 
land. 



[1] H. Fritzsch, Phys. Lett. B 224, 423 (1989). 

[2] C. Amsler, M. Doser, M. Antonelh, D. Asner, K. Babu, 

H. Baer, H. Band, et al., Phys. Lett. B 667, 1 (2008). 
[3] J. Aguilar-Saavedra, Acta Phys. Polon. B35, 2695 

(2004), hep-ph/0409342. 
[4] F. Larios, R. Martinez, and M. A. Perez, Phys. Rev. D 

72, 057504 (2005), hep-ph/0412222. 
[5] P. Fox, Z. Ligeti, M. Papucci, G. Perez, and M. D. 

Schwartz, p. 2 (2007), arXiv:0704.1482 [hep-ex]. 
[6] The CDF-II detector is described in more detail in many 

recent publications; see, for example, A. Abulencia et al., 

(CDF Collaboration), Phys. Rev. D 73, 112006 (2006), 

and references therein. 
[7] N. Cabbibo, Third Topical Conference on Proton- 

Antiproton Collider Physics, Rome, Jan. 1983, CERN 

83-04, Geneva, p. 567 (983); F. Halzen and M. Marsula, 

Phys. Rev. Lett. 51, 857 (1983); K. Hikasa, Phys. Rev. D 

29, 1939 (1984); N.G. Deshpande et al., Phys. Rev. Lett. 

54, 1757 (1985); A.D. Martin, R.G. Roberts, and W. J. 

Stirling, Phys. Lett. 189B, 220 (1987): E.L. Berger, F. 

Halzen, C.S. Kim, and S. Willenbrock, Phys. Rev. D 40, 

83 (1989). 

[8] E. Abouzaid and H. Frisch, Phys. Rev. D 68, 033014 
(2003). 

[9] B. Cooper, Ph.D. thesis. University CoU. London (2006), 

FERMILAB-THESIS-2006-61. 
[10] A. Messina (CDF Collaboration), Braz. J. Phys. 37, 840 

(2007), arXiv:0708.1380 [hep-ex]. 
[11] Transverse momentum and energy are defined as px = 

psin^ and Et = iJsin^, respectively. Missing Et (-^t) 

is defined by = — E^hi, where i is the calorimeter 



tower number for \rj\ < 3.6 (see Ref. [TH]), and hi is a unit 
vector perpendicular to the beam axis and pointing at the 
i**" tower. We correct for jets, photons, electrons and 
muons. We define the magnitude JErp = \^t\- We use the 
convention that "momentum" refers to pc and "mass" to 
mc^. 

[12] F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 80, 
2525 (1998). 

[13] LEP Electroweak Working Group (EWWG), 
|http : //lepewwg . web . cern . ch/LEPEWWG/ [ 

[14] T. Aaltonen et al. (CDF Colaboration), Phys. Rev. Lett. 
101, 192002 (2008), arXiv:0805.2109 [hep-ph]. 

[15] The CDF coordinate system is cylindrical and right- 
handed, with the X axis horizontal and out of the Teva- 
tron ring, t he y axis up, the z axis along the proton beam, 
and r = y'x^+j/^. ^ jg ^j^g azimuthal angle. The pseu- 
dorapidity 77 = — ln(tan(S/2)), where 9 is the polar angle; 
the rj regions for detector components are defined with 
respect to the center of the detector. The central region 
is defined as \rj\ < 1. 

[16] A. Sill et al., Nucl. lustrum. Meth. A 447, 1 (2000); A. 
Affolder et al., Nucl. lustrum. Meth. A 453, 84 (2000); 
C.S. Hill, Nucl. lustrum. Meth. A 530, 1 (2000). 

[17] A. Affolder et al., Nucl. Instrum. Meth. A 526, 249 
(2004). 

[18] T. Aahonen et al. (CDF Collaboration), Phys. Rev. D 

77, 112001 (pages 48) (2008). 
[19] L. Balka et al., Nucl. Instrum. Meth. A 267, 272 (1988). 
[20] S. Kuhlmann et al., Nucl. Instrum. Meth. A 518, 39 

(2004). 

[21] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71, 



24 



032001 (2005). 

[22] M. G. Albrow et al. (CDF Collaboration), Nucl. Instrum. 
Meth. A 480, 524 (2002). 

[23] F. Abe et al., Phys. Rev. Lett. 68, 1104 (1992). 

[24] G. Ascoli et al. (CDF Collaboration), Nucl. Instrum. 
Meth. A 268, 33 (1988). 

[25] T. Dorigo et al. (CDF Collaboration), Nucl. Instrum. 
Meth. A 461, 560 (2001). 

[26] D. Acosta et al. (CDF Collaboration), Nucl. Instrum. 
Meth. A 494, 57 (2002). 

[27] F. A. Berends, W. T. Giele, H. Kuijf, R. Kleiss, and W. J. 
Stirling, Phys. Lett. B 224, 237 (1989). 

[28] P. Renton (2008), arXiv:0804.4779 [hep-phj. 

[29] A. Abulencia et al. (CDF Collaboration), J. Phys. G34, 
2457 (2007), hep-ex/0508029. 

[30] For central electrons at least 5 hits in each of 3 axial and 
2 stereo layers of the COT are required. 

[31] The fraction of electromagnetic energy allowed to leak 
into the hadron compartment Ehad/Ecm must be less 
than 0.055+0.00045 x Ecm (GeV) for central electrons, less 
than 0.05 for electrons in the end-plug calorimeters, less 
than max[0.125, 0.055+0.00045 x Ecm (GeV)] for photons. 

[32] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71, 
051104 (2005); hep-ex/0501023. 

[33] For tight muons at least 5 hits in each of 3 axial and 3 
stereo layers of the COT are required; for loose muons 
with a matching muon stub this is relaxed to 3 axial and 
2 stereo. Loose muons without a matching stub have an 
additional cut on the for Ihe fit to the track. 

[34] The energy deposited in the calorimeter tower traversed 
by the muon must be less than 2 + max(0, 0.0115 x (p — 
100)) GeV in the electromagnetic compartment and less 
than 6 + max(0, 0.028 x (p - 100)) GeV in the hadronic 
compartment. 

[35] A. V. Varganov, Ph.D. thesis. University of Michigan 
(2004), FERMILAB-THESIS-2004-39. 

[36] The muon "stub" in the muon systems must be within 7, 
5, and 6 cm of the extrapolated COT track position, in 
the CMU, CMP, and CMX muon systems, respectively. 

[37] A. V. Kotwal, H. K. Gerberich, and C.Hays, Nucl. In- 
strum. Meth. A 506, 110 (2003), ISSN 0168-9002. 

[38] The maximum correction to the identification efficiency 
for central electrons is 1.7%, for "plug" electrons 6.3%, 
and for central muons is 7.4%. 

[39] F. Abe et al. (CDF Collaboration), Phys. Rev. D 45, 
1448 (1992). 

[40] A. Bhatti et al., Nucl. Instrum. Meth. A 566, 375 (2006), 

hep-ex/0510047. 
[41] A. A. Affolder et al. (CDF Collaboration), Phys. Rev. D 

64, 032002 (2001), hep-ex/0101036. 
[42] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71, 

052003 (2005). 

[43] The transverse mass is defined as Mtrans(^i') ~ 
y/2pT{l)pT{i')il - cos {(j){£) - ct>{v))), where Pt(X) is 
the transverse momentum of X {X is £ or v) and 4>{X) 
is the azimuthal angle of X. 

[44] R. Field (for the CDF Collaboration), AIP Conf. Proc. 



828, 163 (2006), we use Tune A. 
[45] We use Version 2.10 of alpgen. 

[46] PYTHIA v6.216 with a Mi'y'/Z) > 30 GeV cut. Tune A, 
and the "Willis Sakumoto" corrections applied [63|. The 
Z + HF samples are produced with ALPGEN v2.10-prime 
(which has built-in MLM matching [6l[) and showered 
with PYTHIA V 6.326. The both generators used CTEQ5L 
parametrization of the parton distribution functions. 

[47] J. M. Campbell and R. K. Ellis, Phys. Rev. D 60, 113006 
(1999). 

[48] The theoretical estimate of the top pair production cross- 
section is based on (5^ using a top-quark mass of 170.9 
± 1.8 GeV m. 

[49] J. Alwall et al., Eur. Phys. J. C53, 473 (2008), 
arXiv:0706.2569 [hep-ph]. 

[50] The scale factor for W+HF processes (W + bb, W + cc, 
and W + c) does not depend strongly on the number of 
ti events since the tt — > WbZc and tt — > WbWb modes 
contribute mostly to a final state with a W and three or 
four jets in the final state. 

[51] The fit for W — > fj,v + 4 jets with at least one jet b-tagged 
returns zero events with a large uncertainty, due to the 
small statistics of the 4-jet sample. We consequently use 
the value of 2.6% obtained from the fit for inclusive W —> 
fiv + 4 jets, consistent with the fit. 

[52] J. AlwaU et al., JHEP 09, 028 (2007), arXiv:0706.2334 
[hep-ph]. 

[53] A. Abulencia et al. (CDF Collaboration), Phys. Rev. D 

73, 032003 (2006), hep-ex/0510048. 

[54] A. Abulencia et al. (CDF Collaboration), Phys. Rev. D 

74, 072006 (2006), hep-ex/0607035. 

[55] M. L. Mangano, M. Moretti, F. Piccinini, R. Pittau, and 
A. D. Polosa, JHEP 0307, 001 (2003), hep-ph/0206293. 

[56] T. Sjostrand, L. Lonnblad, and S. Mrenna (2001), hep- 
ph/0108264. 

[57] The complete formulas for the expected numbers of 
events are: 

Xi^ = Bi^ + Ntt ■ {(1 — Bz)^ ■ Aww^i + 
Bz{l — Bz) ■ {Awz^e g^) + 
bI-Azz^iVt} (31) 

and 

Xu = Bee + Nu ■ {Azw^te ■ Bz + 

{Azz^u — Azw^u) ■ B'z}. (32) 

[58] M. Cacciari, S. Frixione, G. Ridolfi, M. L. Mangano, and 

P. Nason, JHEP 04, 068 (2004), hep-ph/0303085. 
[59] CDF and DO Collaborations (2007), Fermilab-TM-2380- 

E, TEVEWWG/top 2007/01, hep-ex/0703034. 
[60] D. Acosta et al. (CDF Collaboration), Phys. Rev. Lett. 

94, 091803 (2005), hep-ex/0406078. 
[61] M. L. Mangano, M. Moretti, F. Piccinini, and M. Trec- 

cani, JHEP 01, 013 (2007), hep-ph/0611129. 



